
/* 
Placebo regressions
*/

************
*** OLS ****
************

set more off

estimates clear
clear

use ".\Dropbox\Vaccination\covid.dta" // load the data

drop if statefips==48|statefips==15 // drop states without the vaccination data

zscore(flu rep*2016 vacc*)  // normalize the variables

label variable z_rep_pres2016 "Trump vote share, 2016"

local county "lat lon temp rain income native_share white_share col_share elder_share male_share  mfg_share_diff popden" // county characteristics


// without county controls
reg z_flu z_rep_pres2016  i.statefips, cluster(statefips) 
estimates store ols1

reg z_flu z_rep_pres2016 `county' i.statefips, cluster(statefips) 
estimates store ols2

reg z_vaccinated_18pluspop_ratio z_rep_pres2016  i.statefips, cluster(statefips) 
estimates store ols3

reg z_vaccinated_18pluspop_ratio z_rep_pres2016 `county' i.statefips, cluster(statefips) 
estimates store ols4

estout ols* using ".\Dropbox\Vaccination\Draft\tab_placebo_ols.tex",  replace style(tex) cells(b(fmt(3) star) se(par fmt(3))) stats(r2 N, fmt(2 0) labels("R-squared" "Observations")) keep(z_rep_pres2016) label mlabels(none) collabels(none) starlevels(* 0.10 ** 0.05 *** 0.01)
